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ABSTRACT 

Ouyed et al. (1998) proposed Deuterium (DD) fusion at the core-mantle interface of giant planets as 
a mechanism to explain their observed heat excess. But rather high interior temperatures (~ 10® K) 
and a stratified D layer are needed, making such a scenario unlikely. In this paper, we re-examine DD 
fusion, with the addition of screening effects pertinent to a deuterated core containing ice and some 
heavy elements. This alleviates the extreme temperature constraint and removes the requirement of a 
stratified D layer. As an application, we propose that, if their core temperatures are a few times 10^ K 
and core composition is chemically inhomogeneous, the observed inflated size of some giant exoplanets 
(“hot Jupiters”) may be linked to screened DD fusion occurring deep in the interior. Application of 
an analytic evolution model suggests that the amount of inflation from this effect can be important if 
there is sufficient rock-ice in the core, making DD fusion an effective extra internal energy source for 
radius inflation. The mechanism of screened DD fusion, operating in the above temperature range, 
is generally consistent with the trend in radius anomaly with planetary equilibrium temperature Tgq, 
and also depends on planetary mass. Although we do not consider the effect of incident stellar flux, 
we expect that a minimum level of irradiation is necessary to trigger core erosion and subsequent DD 
fusion inside the planet. Since DD fusion is quite sensitive to the screening potential inferred from 
laboratory experiments, observations of inflated hot Jupiters may help constrain screening effects in 
the cores of giant planets. 

Subject headings: planetary systems - planets and satellites: general - planets and satellites: interiors 


1. INTRODUCTION 

The natural assumption in evolutionary models of 
Jupiter-size exoplanets that are as old as their stellar 
hosts (several Gyr) is that they should have evolved to 
their degeneracy dominated configuration, their radius 
converging to a Jupiter-like radius Rj (e.g. Zapolsky & 
Salpeter 1969; Stevenson 1982; Hubbard 1984). But the 
observed size of such exoplanets belies this modeling as¬ 
sumption, even when strong irradiation from the parent 
star is taken into account^. This “radius anomaly” of 
hot Jupiters (Guillot et al. 2006), combined with the 
variation in radii of most known exoplanets, may be con¬ 
sidered a challenge to standard evolutionary models of 
these objects. 

A few explanations put forward are : (i) capture of a 
small percentage of the stellar radiation, or wind, that 
is deposited deep inside the hot exoplanet, providing an 
“internal heat” source to slow down the cooling suffi¬ 
ciently to inflate the planet; (e.g., Baraffe et al. 2010); 
(ii) tidal heating due to unseen companions pumping up 
the eccentricity (Bodenheimer et al. 2001; see also Jack- 
son et al. 2008; Hansen 2008; Miller et al. 2009; Ibgui 
& Burrows 2009); (hi) kinetic heating due to the break¬ 
ing of atmospheric waves (Guillot & Showman 2002) or 
magneto-viscous heating from ionized winds in the mag¬ 
netized upper atmosphere (Batygin & Stevenson 2010); 
(iv) enhanced atmospheric opacity (Burrows et al. 2007); 

^ For Jupiter-mass planets parked at 0.02 AU, radii of ~ 1.2Rj 
are expected at Gyr ages (e.g Burrows et al. 2007). Nevertheless, 
many planets have radii in excess of 1.2Rj, and most receive irra¬ 
diation below that expected at 0.02 AU (e.g. Miller et al. 2009). 


(v) layered convection induced by compositional inhomo¬ 
geneity which slows the contraction due to inefficient heat 
transport in the interior (Chabrier & Baraffe 2007). A 
review of these ideas can be found in Fortney & Nettel- 
man (2010). Huang & Gumming (2012) found it chal¬ 
lenging to inflate massive planets or those with strong 
irradiation using Ohmic heating alone. Burial of heat by 
turbulence (Youdin & Mitchell 2010), and strong winds 
(Ghabrier et al. 2011) are other mechanisms invoked to 
diffuse stellar irradiation into the planet interior, though 
recently, Kurokawa & Inutsuka (2015) argued that lay¬ 
ered convection with a monotonic chemical composition 
gradient may be insufficient to explain the amount of ra¬ 
dius inflation. It has also been proposed that very short 
period, low-mass, binaries could be the progenitors of 
stellar mergers which may lead to the observed popula¬ 
tion of very hot Jupiters (Martin et al. 2011). In this 
model, inflation is an indication of youth. 

As a result, while the radius anomaly of hot jupiters 
probably originates from additional sources of energy dis¬ 
sipation, the details of the mechanism for radius inflation 
are a topic of current debate. 

In light of this, we are motivated to explore in more 
detail the Deuterium fusion model of Ouyed et al. (1998; 
hereafter OFCS98), which we present here as a feasible 
alternative scenario that may be an important factor in 
(or even at the root of) the radius anomaly. As this is a 
truly deep internal nuclear source of energy generation, 
not one that traps the external irradiation through at¬ 
mospheric phenomena, our idea is novel within the range 
of present explanations for inflated hot jupiters. 

In OFCS98, we pursued the possibility that the excess 
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heat of Jupiter might be due to DD fusion in a sedi¬ 
mented layer of pure D only a few km in thickness and 
concentration ^ 2-4 g cm“^. Extrapolating the DD fu¬ 
sion rates from laboratory energies down to the scale of 
planetary cores, we found that a temperature as high as 
T ^ 10® K is needed to explain the excess heat, which 
is an order of magnitude larger than normally expected 
core temperatures. Furthermore, the sedimentation and 
stratification of the D layer required an unspecified mech¬ 
anism. The stability of the layer to convection also posed 
a challenge to this picture. Here, we present a modifica¬ 
tion to the original idea in OFCS98 that alleviates these 
issues and apply it to Jovian exoplanets. By shifting the 
source of D from the inner layers of the envelope to the 
deuterated core and including an effective screening po¬ 
tential, we argue that energy release from DD fusion at 
typical core temperatures of ~ 10^ K is still sufficient to 
slow down the cooling of the exoplanet. With resulting 
luminosities as large as Pdd ~ 10^® erg s“^, screened 
DD fusion in the deuterated core can provide a viable 
persistent, energy source to inflate close-in, hot Jupiters. 

Let us recall that a viable mechanism to explain the 
excess heat and radius anomaly should be consistent with 
the following key facts: 

(i) The energy release must be persistent and must 
be operating on timescales much larger than the Kelvin- 
Helmholtz times (tkh ~ 10® yrs; the time it takes a 
Jupiter-like planet to settle back to its degeneracy dom¬ 
inated configuration); 

(ii) The data shows that the radius anomaly 
AP=i?obs-Ppred (the difference between the observed 
and predicted radius) scales with the planetary equilib¬ 
rium temperature Tgq, with exponents that vary with the 
mass range. While Laughlin et al. (2011) initially found 
a global best-fit dependence Ai? oc Tg“ witha=1.4± 0.6, 
subsequent work (Enoch et al. 2012; Weiss et al. (2013)) 
has found that taking planet mass, incident flux (or or¬ 
bital separation) into account provides more accurate fits 
for subsets of the data. Weiss et al. (2013) have argued 
that incident flux is a better indicator of radius than 
planetary mass for jupiter-mass planets, and their fits 
support Ai? oc Tg“ with Qf=0.5. Their data also show 
that for smaller mass planets, the radius does depend 
strongly on the mass. A proposed mechanism should in¬ 
clude and investigate this mass-radius relationship. 
Guillot et al (2006) had proposed that the radius- 
inflation mechanism is anti-correlated with the parent 
star metallicity, i.e, the radius anomaly grows smaller 
with increasing metallicity. However, this is not a strong 
constraint since if the planets around metal-rich stars 
are themselves metal-rich, they would be naturally be 
denser and more compact (Johnson et al. (2010); Enoch 
et al. (2012)). The issue is still debated as the correla¬ 
tion between parent star metallicity and planet metallic¬ 
ity is weaker than originally thought (Thorngren et al. 
(2015)). We therefore do not take the metallicity trends 
as a significant constraint on the inflation mechanism. 

The above two points can be obtained in our model if 
we restrict the parameter space of the screening potential 
based on the astrophysics. This paper is organized as fol¬ 
lows: In §2 we describe screened DD burning under con¬ 
ditions prevailing in the interior of hot, close-in, Jupiters 
and estimate the screened reaction rate at relevant tem¬ 


peratures. In §3, we apply the results to the core of hot 
jupiters, estimating the fusion power from this reaction, 
the amount of radius inflation from a simple evolution 
model, timescale of the effect, the trend with planetary 
equilibrium temperature Teq and the trend with plane¬ 
tary mass. Concluding remarks are presented in §4. 

2. DD FUSION IN JOVIAN EXOPLANETS 
2.1. DD fusion in a D layer 

In OFCS98, assuming a concentrated layer of D near 
or above the core in giant planets, we showed that DD 
burning can be triggered and sustained if sufficiently high 
temperatures {T 10® K) are maintained once the D 
layer forms^. At T < 10® K (i.e. T 6 eV), the 
DD cross-section decreases sharply (Fig. 1 in OFCS98). 
These are extreme conditions since (i) a 10® K tempera¬ 
ture may barely be achieved, if at all, only early in the 
planet’s history; (ii) the formation, stability and preser¬ 
vation of the D layer is compromised by the strong con¬ 
vective mixing in the interior. In OFCS98, we therefore 
envisaged the D layer to be situated in a non-convective 
portion of the planet’s interior. At a temperature T, 
the reactivity (av) used in OFCS98 is the Bosch-Hale fit 
(Bosch & Hale 1992) based on the functional form sug¬ 
gested by Peres (1979): {av) = cm® s“^ 

where 9 and f are given by eqs. (6) and (7) respectively 
in DECS 98® while yi is the reduced mass of the reactants 
and c the speed of light. The exponential fall of the re¬ 
activity, some 10 orders of magnitude for a factor of 2 
in temperature, is moderated by standard electron and 
ion screening effects so that the effective reaction rate 
Ceff = fefiicrv). Here, fe{T,p) is the strong screening 
limit of the degenerate quantum electron gas while /i(T) 
accounts for screening by the classical ionic liquid formed 
by D nuclei (eqs. (8) & (9) respectively in OFCS98). The 
power generated from the D layer was found to be of the 
order 10^^ erg s“^ for temperatures of ^ 10® K at den¬ 
sities prevailing in the core-mantle region inside planets 
such as Jupiter. Once DD burning starts, the reactivity 
and thus the temperature in the D layer increases. Two 
possible outcomes are: (i) a runaway burning because 
of the degenerate core conditions. In this case, most of 
the D is consumed on very short timescales (compared 
to the planet’s lifetime) leading to a “sudden” release 
of energy; (ii) a steady-state burning at an equilibrium 
temperature (Tc ~ 10® K) below the convection thresh¬ 
old. This second scenario may be a viable explanation of 
the heat excess, but the high temperature and survival 
of the D layer once burning starts are unfavorable and 
adhoc conditions (Ouyed 2004). 

2.2. DD fusion in the deuterated core 

To find a more plausible alternative to these extreme 
conditions, we take a closer look at DD screening in a 
deuterated core environment. Here, DD burning is oc¬ 
curring in the core which is bombarded by the D ions 

^ The Deuterium could be supplied by planetesimal vaporization 
(Ouyed 2004) or, as we argue in this paper, by core erosion (e.g. 
Guillot et al. 2004). 

® Eq.(6) in OFCS98 should read 9 = - T(C 2 +T(Ci+TCe )) — 

1 + T(C3+T(C5+TC7)) 

with T in keV and fitting parameters Ci from Bosch Hale (1992). 
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freed from the icy D-bearing core material by erosion 
(e.g. Guillot et al. 2004). The freed D which is lifted 
upward by convection eventually settles back down onto 
the core triggering screened fusion. Thus we rely on the 
presence of icy material in the core to supply the free D 
(the projectiles or beam) which fuse with the D in the 
core (the target). Other sources and supply mechanisms 
of free D (which we do not consider here) originate in 
the overlaying H-rich envelope and deposition of D deep 
inside the planet’s interior by planetesimal impact and 
vaporization (Ouyed 2004). 

Let us represent the screened cross-section by dscr- 
Several laboratory-based studies have shown that the DD 
fusion cross-section is enhanced in a solid deuterated tar¬ 
get as compared to a gas target (Czerski et al. 2001, 
Raiola et al. 2002, Kasagi et al. 2002, Kasagi 2004). 
This is attributed to correlations between the conduc¬ 
tion electrons in a metal lattice and enhanced mobility 
of deuterons. The environment of the target nucleus, 
in particular, the mobility of D ions in the metal lattice, 
plays an important role in fusion well below the Coulomb 
barrier, and is the principal reason for the large enhance¬ 
ment. The enhancement factor at energy E is (Raiola et 
al. 2002, Kasagi et al. 2004) 


as 

POO 

(ascrv) = / ascT{v)vf{v)dv , ( 2 ) 

JQ 

where v = ■^2E/^ and f{v) is the Maxwellian distri¬ 
bution of speeds. The integral is usually carried out 
in the energy domain by approximating the integrand 
as a narrow Gaussian around the Gamow peak. The 
S-factor is roughly energy independent, and the ther¬ 
mally averaged rate can then be extrapolated to tem¬ 
peratures of astrophysical interest based on fits to multi¬ 
channel cross-sections from R-matrix theory (Bosch & 
Hale 1992). In the case of extrasolar Jupiters, core tem¬ 
peratures are in the range of a few eVs, much lower than 
in stellar cores where temperatures are of order keV. 
A naive extrapolation of the Bosch-Hale fit to fee’T ~ 
eV shows that 9 <T (see footnote 3), so the rate falls 
off at least as rapidly as T“^/^exp(—diHo/fcsT) where 
Eo={bkBT/2y/^; here fee is the Boltzmann constant. 
However, now the screening enhancement boosts the re¬ 
action rate by several orders of magnitude at low tem¬ 
perature. This effect is immediately apparent when we 
approximate E/{E U) ~ Ep\^/U in the pre-factor of 
Eq. O where Epk is obtained from 


CTscr _ E Exp(6£' - 1) 

CTbare ~ E+U Exp(6(£; + - 1) ’ ^ ’ 


for S{E) Ri S{E + U) where b = TTae\/2fic^ with Oe the 
fine structure constant, U is the screening potential for 
the metal target, and S the astrophysical S-factor. The 
value of the screening potential U has a large impact 
on our calculations, since the DD fusion cross-section 
is exponentially sensitive to it at low temperatures. 
Laboratory values for U vary widely from one target 
to the next, even within the same experimental setup 
(eg., from < 30 eV for Al target to ~ 450 ± 50 eV for 
Fe target; Raiola et al. 2002) and have not achieved 
consensus between different groups (Kasagi 2004 reports 
U ~ 200 ± 20 eV for Fe), so that we cannot choose a 
value with confidence a priori. In addition, the absence 
of detailed knowledge of the core composition affects the 
choice of U. While laboratory experiments do extremely 
well to achieve a few hundred eV precision for nuclear 
physics (MeV scales), this is insufficient to constrain 
our model. We will therefore allow the astrophysics 
to further constrain the range of U by requiring that 
the resulting power from DD fusion Pdd allows some 
amount of inflation while lasting for Giga-years. We 
emphasize that the sensitivity to U is an opportunity 
to use exoplanet observations to constrain the screening 
effects in a unique regime of temperature and density. 
Adopting a typical value of U of 200(±20eV) below 
makes our model consistent with observational numbers 
and trends (see sec.§3), and this value corresponds to 
the screening parameter in a deuterated Fe target, as 
per Kasagi et al. (2004). Since D is expected to be 
found in ice pockets, this could be interpreted to mean 
that Fe is the main component of the rocky material 
that diffuses from deeper layers and gets mixed with ice. 

Taking screening into account, we can write the corre¬ 
sponding thermally averaged cross-section of DD fusion 


b ^ E 
y/E + U ^ 

We find that Epk = Eq — U, which adds a new mul¬ 
tiplicative factor exp(U/T) in the Gaussian approxima¬ 
tion to the rate - this explains the strong sensitivity of 
the cross-section to the value of U apparent from Fig. [TJ 
Since ifpk > 0 implies Eq > U, our estimate of the en¬ 
hancement factor is accurate for temperatures that are 
large enough to satisfy fc^T > U^/‘^{2/b) (typically 5 eV 
or more if we choose U ~ 200 eV). Taking {7 —>■ 0 yields 
the usual temperature dependence of the reactivity in the 
absence of screening: exp(—OEo/fesT). But core temper¬ 
atures in extrasolar Jupiters can be somewhat lower than 
5eV, so we make a more accurate calculation of Epk valid 
for Eq <U hy setting 

d f b In 

The Gamow peak is now located quite close to ksT, 
which tempers the exponential fall-off with temperature. 
As before, we perform a Gaussian approximation to the 
integrand to evaluate the integral. The reaction rate 

{ctsctv) now behaves as ^ r^/^exp(—6/-\/c7)xexp 

While this is a monotonically increasing function of 
temperature, the effective cross-section aeS=fi{o'sciv) is 
monotonically decreasing because of the strong effect of 
ionic screening fi{T) at low temperatures. The factor 
fe{p,T) need not be included in since the screen¬ 
ing due to conduction electrons is already present in 
U. The two approximate expressions for the rate tTeff, 
one rising and the other falling with T on either side 
of keT r; U^/^{2/b), are joined together by a cubic 
Hermite spline to obtain a smooth continuous curve for 
the rate around a few eVs. The two individual rates 
are naturally quite sensitive to C/, therefore, we expect 
the fitted fusion reactivity to depend strongly on U as 



E + U 


E 


= 0 . (4) 



= 0 . 


( 3 ) 




pk 
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Fig. 1.— Reactivity (o-scrii) (solid lines) from Eq. l[2} for 
different values of screening potential U from 100-300 eV, and the 
corresponding effective cross-section (Teff = fl{T)(crscTv) (dotted 
line) for 5 X 10® K < T < 5 X 10'* K, representative of the core of 
giant planets. The value of 17=200±20 eV for Eq. |[TJ corresponds 
to an Fe target (Kasagi 2004) and appears to fit best with our 
model. 


Fig. 2.— Fusion power Log(PQD) from Eq. ([TJl for different val¬ 
ues of screening potential U (in eV) for 5 x 10^ K < T < 5 x 10^ K, 
representative of the core of giant planets. A value oi U ^ 180- 
220 eV is appropriate for Fe-dominated rock-ice, which gives 
^DD ~ lO^^-lO^® erg s~^, yielding sufficient inflation. 


well, and indeed this is what is seen in Fig. [H which 
shows the reactivity {cscrv) with the effective potential 
U inclnded, alongside the corresponding effective cross- 
section deff = fiiT) (dscr d) for 5 X10^ K < T < 5 X10^ K, 
representative of the core of giant planets. 

3. APPLICATION TO GIANT/HOT PLANETS/JUPITERS 
3.1. Fusion Power 

Taking into account the two main channels of DD fu¬ 
sion, viz., d(d,n)^He and d(d,p)t (see the DD cycle and 
Fig.3 in OFCS (1998)), we compute the power release 
Pdd for a specific core density of D (pd=ii-d?iid) and U 
value"^ as follows (see Eq. (14) in OFCS98): 

DDD:i(C^) = [ Qi^CTeff(C/)47rr^dr- , (5) 

Jo ^ 


core-mantle interface due to sedimentation or core ero¬ 
sion (Wilson & Militzer 2012a; 2012b) leading to chemi¬ 
cally inhomogeneous layers deep in the planetary interior. 
In these works, it was not clear that the solvated mate¬ 
rials could be convected upwards efficiently through the 
lighter layers on top. However, Wilson (2015) has pro¬ 
posed that the mechanism of semi-convection can still 
lead to high diffusivity of the heavy rock-ice into the up¬ 
per layers, based on modern density-functional calcula¬ 
tions. Semi-convection separates layers of constant com¬ 
position with thin diffusive layers that are well-mixed, 
leading to a slower rate of redistribution and a slower 
rate of cooling. This also implies that core erosion and 
mixing can persist for longer times. In such a scenario, 
our model becomes viable since DD fusion will occur in a 
deuterated metallic environment, supported by slow core 
erosion and slow cooling. Assuming a uniform spatial 
profile for the core density, we have 


where is the core radius and wd is deuterium mass. 
The index i = 1,2 stands for the two main channels of 
DD fusion with Qi the energy released per D atom fused 
for each channel. In terms of core density pc, the core 
density of D is pD = 2(D/i7)xpH = 2(D/i7) x (2/18)pice, 
which gives 

PD -- 2.2 X 10“®?7d -4?7ice.-lPc , (6) 

where pice,-i is the amount in mass of ice in the core 
in units of 10“^ (i.e. pice=0.1pc) and pd ,-4 the D to H 
ratio (D/H) in the core in units of 10“"^. The value of 
D/H ratio chosen is a few times higher than its protoso¬ 
lar value, since this isotopic ratio is higher in planetary 
ice than interstellar medium (Guillot (1999)). Within 
our model, one must address the likelihood that heavy 
elements from the rocky, icy core can be well-mixed with 
the less dense layers of H above it, so that the D liberated 
by core erosion finds itself in a metallic deuterated envi¬ 
ronment. Recent works have argued that heavy elements 
from the core can be mixed with fluid H/He layers at the 

Although Pdd depends weakly on temperature T through the 
reactivity cTeffi we do not write it explicitly since the dependence 
with U is more important in the temperature range 10^-10®K. In 
fact, Pdd for our median value oi U ~ 200 eV is almost constant 
in this temperature range (see Fig. 0. 


PDD;i(C/) ~ ^(2.2xl0-6)2,y2 4^2 -^i^aeff(C/) , 

2 ’ ’ rriY) rriY) 

( 7 ) 

where Me = /q ° pedtr-r^dr is the core mass. 

The amount of energy released is Qdd = Q 1 +Q 2 — 2x 
10“® ergs per D atom fused, with Pdd = Pdd,i + Pdd., 2 - 
The mass of the core is hereafter taken to be 10% of the 
planet’s mass (Mc=0.1Mp). For an evolved Jupiter size 
(Pj) and Jupiter mass (Mj) planet, with pc,j ~ 4.38 
g cm“^ as core density (e.g. Hubbard 1984), we find 
Pdd(C/) (3.8 X 1026^) The 

effective fusion reactivity aes{U) is scaled in units of its 
typical magnitude 10“^^ cm^ s“^ for f/=180 eV (typical 
of rock-ice dominated by Fe), and a burning temperature 
in the ~ 10^ K range. We can get an idea of the 
sensitivity of Pdd{U) to the screening potential in Fig.[^ 
which shows the result for Pdd as a function of core 
temperature. 

3.2. Amount of Radius Inflation 

To study the amount of inflation that results from DD 
fusion over Gyr timescales, it is useful to employ a one¬ 
dimensional analytic time-evolution model for the plan¬ 
etary radius Pp(t) (e.g. Guillot 2005). I.e. dE/dt = 
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Fig. 3.— Time evolution of Rp{t) including DD fusion in the 
planetary core, assuming rjice = 0.1 and no stellar irradiation. 
Masses Mp are in units of Jupiter mass. Sub-jupiter mass planets 
are inflated the most, due to smaller heat loss (Tg^; see ^3.51 1 
while higher mass planets have larger internal heat from Pdd fo 
balance cooling. 

Pooit) — Leg{Tee) where E = —3GMp/7Rp{t) is the 
(time-dependent) internal energy (we make use of the 
virial theorem to link the planet’s internal energy to grav¬ 
itational energy), and Leff('7eff) = 47ri?p(t)^crsB'7)fff is the 
standard blackbody luminosity; ctsb being the Stefan- 
Boltzmann constant. The time evolution of Puuit), 
which is the additional internal heat source, comes from 
the following factors: the core mass Me and density pc, 
and the effective cross section OeS which depends on the 
core temperature Tc. We note from Fig.[T]that although 
the core temperature may change by a factor of a few 
during the contraction, the cross-section is only weakly 
dependent on temperature for screening potential U typi¬ 
cal of Fe-dominated rock-ice, so that its time-dependence 
can be ignored. 

We express the core density in terms of the planetary 
mass and radius following the prescription given in Sea- 
ger et al. (2007; see Eq.(ll) and Table III in that pa¬ 
per) together with pressure balance between the core 
and the inner envelope. This effectively gives us the 
time-dependence of the core density in terms of dRp/dt. 
The core mass erodes at a rate proportional to the ex¬ 
isting core mass (since erosion powers Pdd) leading to 
Mc(t) = Mc(0)exp(—O.lt) with Mc(0) = O.lMp. As time 
is measured in Gyr, the core mass effectively decreases by 
« 10% on a timescale of 1 Gyr and about 65% in 10 Gyr. 
We use the mean adiabat connecting the isothermal core 
at temperature Tc to the effective Teff which determines 
the cooling flux (ignoring the effect of stellar irradiation 
of the atmosphere). We find that Tes scales as 
for a mean adiabatic index of r=0.3, as commonly used 
for Jupiter interiors (eg. Hubbard et al. 1968). Here, 
Teff = (200 K)(Mp/Mj)0-4. 

The solution of the evolution equation for R{t) is plot¬ 
ted in FiglSl from which we infer that DD fusion can 
provide enough energy to prevent prolonged contraction, 
and ultimately lead to an equilibrium inflated radius that 
persists for longer than the Kelvin-Helmholtz timescale. 
We discuss the trend with planetary mass in 113.51 

Taking into account opacity effects in the atmosphere 
leads to even more inflation. The semi-analytic treat¬ 
ment in Ginzburg & Sari (2015) which expresses AR in 
terms of the deposited luminosity and optical depth for 



Fig. 4. — Scaling of planetary radius (with respect to 
Jupiter radius) with DD fusion power (relative to equilib¬ 
rium luminosity Deq) for different values of the optical depth 
r = {10^, 10®, 10"^, 10®} and corresponding pressure at the base of 
the radiative layer P = {1,10,100, 1000} bars, using Eq. {S)!. 

a deep internal heating source supports our numerical 
estimate from the radius evolution above. If we com¬ 
pare Pdd to Leq=47ri?p(TSBTg^q, ~ 10^®-10^° erg s“^ (the 
planet’s equilibrium luminosity and Tgq. the equilibrium 
temperature), we find that fusion power can be from a 
few percent (for Mp ~ Mj) up to 100% (for Mp ^ lOMj) 
of Leq- This heat, being deposited near the core bound¬ 
ary, can be convected up rapidly through fluid motion 
and deposited at the base of the radiative layer, where 
thermal resistance is large. In this deposition region, 
opacity from negative Hydrogen ion H“, kh, dominates 
the heat transfer at Teq of a few thousand Kelvins (e.g. 
Arras & Bildsten 2006). Using Eq. (31) in Ginzburg & 
Sari (2015), we hnd 

AR - 0.3i?j (^1 + , (8) 

where the factor 0.3 accounts for inflation from stellar 
radiation alone and a ^ 0.2. This form is in good agree¬ 
ment with the more detailed numerics in Spiegel & Bur¬ 
rows (2013) for arbitrary internal heat sources located 
deep inside the planet. The optical depth r ^ unP/g 
(with g = GMp/Rp) can reach values as high as 10^- 
10® for Teq. > 2000 K in a regime where the H“ opac¬ 
ity dominates. Accordingly, FigH] shows that our model 
provides conditions (i.e. the Pdd/L eq ratio) that allow 
variable amounts of inflation, upwards of 50% and even 
as much as 100%, depending on the optical depth. We 
have also verified that this estimate is in good agreement 
with Figure 6 in Miller et al. (2009), which is also based 
on an internal heating source. Our Fig. |T] is similar to 
Fig. 4 in Ginzburg & Sari (2015) which relies on more 
general arguments, supporting our contention that DD 
fusion can provide enough energy to inflate the planet as 
a true internal heat source. 

3.3. Timescale of Radius Inflation 

As mentioned in the introduction, the inflation should 
persist for times larger than tkh ~ 10® yrs. To esti¬ 
mate the timescale for this effect, we consider the to¬ 
tal energy released by burning a mass Mpt ^ 2.2 x 
10“®?7D,-4Tce,-iAfc of D in the core, corresponding to 
the D density in Eq. (©: 
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iV^D 49 ikfp 

= QdD- ~ 2.6 X 10^^ erg 77D.-4?7ice,-lxr ■ (^) 

Tod -Mj 

This energy will be gradually consumed at a non¬ 
constant rate determined by Pdd- The timescale to 
reach the equilibrium radius Req (where dR/dt « 0) 
gives Pdd ~ Teq = Since we have TeS ^ 

(200K) (Mp/Mj)°-^, we estimate that 

^DD „ p /Mj \ ° ® 

t-dd = -1-0 Gyr ?7 d- 4?7ice-i -jT ’ 

-fDD V-'Wp/ 

where we have used Req=\.2Rj . The above timescale 
of ~ 1 Gyr is an approximate lower bound for Jupiter- 
mass planets, given that the inflation effect can last 
longer as also decreases slowly during the evolu¬ 
tion, though we assumed it constant. This estimate sug¬ 
gests that higher mass planets will have shorter inflation 
timescales, which agrees with FiglS] The timescale also 
scales with ?7ice, which our evolution code also confirms 
- a larger ice content leads to a more persistent inflation 
effect. If ice occurs at the expense of metal content in 
the core, planets with lower metallicity (presumably de¬ 
rived from the parent star) would be more likely to be 
inflated as observed now, having a long-lived source for 
DD fusion. We remark on the observational trend with 
metallicity in (2) We conclude that according to Eq. (ITUl) 
or Fig. [31 Jupiter mass planets would inflate by about 
20-50%, with the effect clearly lasting for times of order 
10® yrs or more, which is larger than Ikh- Besides the 
fusion power, inflation amount and timescale, there are 
other aspects of our model that we would like to explore. 
In the next section, we discuss some interesting trends in 
the observational data that highlight possible sources of 
variability in our model. 

3.4. Dependence of Radius Anomaly on T^q 

Based on observations of over 90 transiting giant ex¬ 
oplanets, Laughlin et al. (2011) showed that the radius 
anomaly for Tgq > 1000 K can be anywhere from 0 to 
70%, and that despite this variability, there appears to be 
a general rising trend between the radius anomaly AR = 
7?obs-7?pred and effective planetary temperature Teq with 
a best-fit dependence AR oc with a = 1.4 ± 0.6. This 
cannot be explained by variations in planetary mass or 
stellar irradiation alone. Laughlin et al. (2011) used this 
correlation to argue for Ohmic dissipation from MHD ef¬ 
fects in surface atmospheric flows as the preferred mech¬ 
anism for inflation over kinetic heating. Since we are 
proposing quite a different mechanism (deep internal 
heating), it is interesting to see whether our model can 
accommodate this scaling. As mentioned in the intro¬ 
duction, subsequent work (Enoch et al. 2012; Bayliss 
et al. (2103); Weiss et al. (2013)) used a larger sam¬ 
ple size and more accurate fits for subsets of the data 
that take other planetary properties such as mass into 
account, leading to different exponents in the scaling re¬ 
lation. Therefore, we also wish to study the dependence 
of the radius inflation on planetary mass. 

Let us first focus on the AR scaling with the plan¬ 
etary equilibrium temperature Teq. The latter charac¬ 
terizes the incoming flux, and (roughly speaking) rep¬ 
resents an average photospheric temperature given by 


Teq = (i?*/(2a*))^^^T* for zero albedo and orbital ec¬ 
centricity, where * denotes stellar numbers and a is the 
semi-major axis (Laughlin et al. 2011). A general way 
to describe radius inflation from a deep energy deposi¬ 
tion, assuming an n = 1 polytrope for Jupiter-like plan¬ 
ets, is given by Eq. JS]), taken from Ginzburg & Sari 
(2015), which agrees well with the numerical estimates 
presented in Spiegel & Burrows (2013). Since the usual 
heat-blanketing effect from stellar irradiation is insuffi¬ 
cient to explain the observed inflation, it follows that 
additional heat trapping must occur. This happens if 
Pddt > Teq and is the effect described by Eq. ([8]). In 
our model, as shown by Fig. 21 we certainly have regimes 
where ToDT/Teq. >> 1. Therefore, Eq. ([5]) becomes 
AR oc (TDDT/Teq.)“. For a Rosseland opacity (i.e. a 
gray approximation) with Teq > 2000 K, the opacity kh 
is dominated by that of the negative Hydrogen ion H~ 
and scales as kh oc with s >> 1 and can be as high 
as s = 9 (e.g. Allard et al. 2001; see also Fig. 1 in Ar¬ 
ras & Bildtsten 2006). Since Teq oc T^q, this means that 

AR oc From Fig. O we see that Tdd is more 

or less independent of Tc, the core temperature. Since 
Tc and Tgff (the temperature at the base of the radiative 
layer) lie on the same adiabat (eg., Hubbard 1977; Guil- 
lot& Showman 2002), we can find a relation between T^ 
and Teq assuming Tgq ~ lOTe® (eg., Spiegel & Burrows 
2013). For an n = 1 polytrope, this yields Tc ~ lOOTgff, 
which in turn implies Tc ^ lOTeq. This linear proportion¬ 
ality implies that Pdd does not scale significantly with 
Teq and that, in our model, the scaling of radius inflation 
with Teq is driven by the optical depth at which the heat 
is deposited, once the condition PoDT/Teq, >> 1 is sat¬ 
isfied. Physically speaking, DD fusion provides a large 
enough power source, deep enough in the interior, to re¬ 
inforce the slowing down of the planet’s cooling due to 
trapping of stellar irradiation. Substituting typical val¬ 
ues for a ^ 0.2 (Ginzburg & Sari (2015)) and taking into 
account the steep dependence of H“ opacity on temper¬ 
ature with s = 9, we get AR cx Tgq®, which agrees within 
systematics to the fit presented in Laughlin et al. (2011). 
We stress that the application of Eq.® to obtain the 
approximate scaling AR cx Tgq^=*=® ® claimed in Laughlin 
et al. (2011) requires TDO^/Teq. >> 1, which in our 
model, holds true for more inflated planets when we ap¬ 
ply the heuristic equations of Ginzburg & Sari (2015). 
However, the fit in Laughlin et al. (2011) works better 
for moderately inflated planets (20% or so). The differ¬ 
ence could be explained by choosing a smaller value than 
s=9 for the H~ opacity, which is on the higher end. A 
smaller opacity leads to less heat trapping and less infla¬ 
tion for the same Tdd and r, while giving a scaling that 
is weaker than AR oc T}„^ but still consistent with the 
bounds in Laughlin et al. (2011) or Weiss et al. (2013). 
We emphasize that the fact that Tdd is at best weakly 
dependent on Teq is important in getting the correlation 
index around 1.0, which further supports the choice of 
U in the right range. Thus, we can be optimistic that 
the astrophysics implies a certain self-consistency in our 
model’s assumptions. 

A related finding which is generally consistent with the 
scaling of AR with Teq is that planets receiving mod¬ 
est irradiation seem to lack inflation (e.g. Demory & 
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Seager 2011). In our model, the condition for core ero¬ 
sion/solubility, Tc > Tero (with Tero ~ 10"^ K; Wilson & 
Militzer (2012a&b)), translates to a limit on the equi¬ 
librium temperature Teq > 10^ K x (Tero/IO^K), or in 
terms of irradiation Fi„=asBTeq > 5.6 x 10^ erg cm“^ 
s“^ X (Tero/IO^^K)^. This constraint on Teq, which is re¬ 
flected in the observations, emerges as a necessary con¬ 
dition in our model but is not sufficient for inflation, 
since planets with higher metallicity (i.e. with reduced 
ice content) and/or smaller core will experience little to 
no inflation from DD fusion. 

3.5. Dependence of Radius Anomaly on Planetary Mass 

Now, let us look at the dependence on planetary mass 
in our model. Data from 119 transiting planets in dif¬ 
ferent mass regions studied in Enoch et al (2012) shows 
that a different scaling of radius with Teq for Saturn- 
mass planets, Jupiter-mass planets and high-mass plan¬ 
ets yielded better fits to the radii of these planets. Bayliss 
et al. (2013) produced a density plot of 125 transit¬ 
ing exoplanets in short-period orbits that revealed an 
under-density of planets in the 1-2-1. 6Mj range. Fig. 
11 of Enoch et al. (2012) and Fig.l of Laughlin et 
al. (2011) suggests higher mass planets are on average 
under-inflated. While more data is needed to clarify the 
trends, planetary mass clearly plays a role. Mass de¬ 
pendence enters into the DD fusion mechanism mainly 
through two factors. The first is the core density which 
depends on planetary mass. We use the fitted core EoS 
from Eq.(ll) in Seager et al. (2007), with fit constants 
appropriate to an Fe-rich core (Table III of Seager et 
al. (2007)). Then, pressure balance relates core pressure 
to the stellar mass and radius which yields the required 
mass dependence of the core density for our evolution 
model. The second factor is the mean adiabat which 
connects the isothermal core to Tgg which determines the 
cooling flux (ignoring the effect of stellar irradiation of 
the atmosphere). For a mean adiabatic index of r=0.3 
(eg. Hubbard et al. 1968) this yields Teff scaling as 
TVfo.4_ account the fact that sub-Jupiter 

mass planets are less compressible than super-Jupiters 
by choosing an appropriate polytropic R-AI relation, i.e, 
n=3l2 for high-mass planets and n=0. 5-1.0 for lighter 
planets. Incorporating this dependence into our evolu¬ 
tion model, FigI3] shows that sub-Jupiter (0.1-0.5 Mj) 
and super-Jupiter mass planets (4-10 Mj) are more in¬ 
flated than planets in the mass range 1-2 Mj. The reason 
for this trend within our model is that smaller mass plan¬ 
ets have smaller Teff which slows down the cooling while 
the larger mass planets have higher core mass and den¬ 
sity leading to a larger Tdd that offsets the larger Teff. 
Comparing to Fig. 11 of Enoch et al. (2012) and Fig. 
1 of Laughlin et al.(2011) suggests higher mass planets 
are on average under-inflated, which is not the case in 
our mechanism. A slight change in the mean adiabat de¬ 
pendence from the fully convective case of n=5/3 does 
lead to under-inflation of large mass planets in our model 
since then Teff scales as or higher. In any case, it 

appears that our model needs to be re-examined for the 
high-mass region. We can conclude that our model cer¬ 
tainly gives some of the right trends with mass, while a 
more comprehensive treatment of the interior structure 
of hot Jupiters needs to be taken up in future work to 


better confront our model with the observational trends. 

4. DISCUSSION & CONCLUSIONS 

In the preceding sections, we have presented a model 
of DD burning in the core of extrasolar Jupiters in an at¬ 
tempt to explain the anomalously inflated radius of these 
planets. In this section, we discuss some implications 
and limitations of the model. Since we have presented 
order of magnitude estimates and used a simple time- 
evolution model for the radius, we do not claim that DD 
fusion is the mechanism for radius inflation. However, 
for hot Jupiters where other mechanisms may prove in¬ 
sufficient to inflate these planets, our estimates clearly 
show that DD fusion in the core could be an additional 
heat source that acts as a “hot plate” inside the simmer¬ 
ing planet, leading to over-inflation of their radius. For 
this work, we assumed that core erosion is induced by 
convection (e.g. Guillot et al. 2004). However, cracking 
and sputtering induced by DD fusion (e.g. Smentkowski 
2000) may also lead to self-sustained core erosion (inde¬ 
pendent of convective plumes eroding the core). It may 
be that convection is needed early on to provide enough 
D-bearing material and then sputtering takes over. Once 
DD is triggered, such a self-regulating process can, in 
principle, take place between erosion (fuel supply) and 
burning (fuel consumption) leading to an equilibrium 
temperature in the 10"^ K range, as required for efficient 
core erosion, according to studies by Wilson & Militzer 
(2012a&b). 

Core erosion is a necessary condition (to supply the 
projectiles) but this may not be sufficient if not enough 
D in the core is available for burning. According to 
Eq. ((61), reduced ice content in the core (i.e. reduced D 
content) leads to lesser heat generation. If the increase 
in core metallicity is at the expense of ice (and thus D 
content), then planets with metal-rich cores should expe¬ 
rience less inflation despite ongoing core erosion. Guil¬ 
lot et al. (2006) did claim such a correlation, but sug¬ 
gested that planets around metal-rich stars are them¬ 
selves metal-rich, making them denser and more compact 
(less inflation). They did not consider the possibility of 
internal heating from an icy core. In our model, metal¬ 
licity can be further correlated with reduced ice content 
(e.g. planets forming inside the snow line) which directly 
impacts the amount of inflation. Combining the sug¬ 
gested explanation of Guillot et al. (2006) with ours, a 
natural connection emerges to the snow line in the disk: 
planets that formed beyond the snow line should carry 
healthy amount of ice in their cores making DD burn¬ 
ing a more likely possibility when they migrate and park 
close to the star, and such hot Jupiters should show the 
most inflation. If the parent star metallicity and planet 
mass is known (eg.. Miller & Fortney (2011), one can 
then hope to use the radius trend with metallicity to 
constrain the planets core metallicity, which is indirectly 
a parameter in our model through rjice- However, the 
issue of metallicity is complex, as even the correlation 
between parent star metallicity and planet metallicity 
has been questioned recently (Thorngren et al.(20I5)), 
while other works Johnson et al. (2010) and Enoch et al. 
(2012) are more circumspect on the trend with metal¬ 
licity than Guillot et al. (2006). Therefore, while our 
model would predict an anti-correlation trend, specially 
if metallicity is primarily residing in the core, we can- 
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Fig. 5. — Time evolution of ijp(f) including DD fusion in the 
planetary core with r]ice = {0.2, O.f, 0.05} and no stellar irradiation 
for a Jupiter-like planet {Mp = Mj). 

not infer any constraints on our model from the current 
trends in the data. 

Inflation seems to scale with the irradiation flux re¬ 
ceived by the planet (e.g. Laughlin et al. 2011), while 
in our model which does not take irradiation into ac¬ 
count, the inflation arises from the high luminosity and 
energy deposition at large optical depth. Since irradia¬ 
tion is not sufficient to explain the observed inflation, the 
observational trend may also be suggesting that hotter 
surface temperature can induce conditions in the core fa¬ 
vorable for erosion to be triggered. In that case, when 
there is sufficient D content in the core, ongoing DD 
burning would help sustain the planet’s inflation and in 
some circumstances over-inflation. In extreme cases (e.g. 
Poo/Leq > 0.01 and r > 10®), inflation can lead to ex¬ 
cessive mass-loss and possibly envelope stripping result¬ 
ing in a naked “super-Earth” core (the eroded core of the 
hot Jupiter). 

The dependence on ice content in our model, implies 
that most inflated planets most likely formed beyond the 
ice line (i.e. accreted the necessary ice along the rocky 
material) and migrated to park close to their parent star. 
Once it collected enough ice in its core, a young (hot) gi¬ 
ant planet is prone to DD burning as long as an erosion 
mechanism is at play. In principle, mechanisms other 
than trapping of stellar radiation (e.g. tidal heating) 
could be capable of inducing conditions in the core that 
trigger core erosion, leading to DD fusion. Thus inflated 
but colder Jupiters (at higher orbits) may exist, but be 
hard to find due to selection effects in transit detections. 
It is thus interesting to speculate on the implications of 
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DD burning in our own Jupiter where the effects of ir- 
radiation can be ignored. Fig. [5] shows the evolution 
of the radius of a Jupiter-like planet as induced by DD 
fusion for different ice content in the core. The conhg- 
uration with less ice content has been slowly shrinking 
to Jupiter’s current configuration purely due to cooling 
and eroding of the core. It suggests that even for planets 
far from their parent star, heat released from possibly 
ongoing DD fusion can keep the core temperature at val¬ 
ues favorable for erosion to proceed and thus sustain the 
fusion. However, these planets will cool faster than the 
hot, parked, Jupiters and should deflate to a radius close 
to the degenerate value after only a few billion years. An¬ 
other possibility, for example, is that planets which, long 
ago, experienced DD fusion while parked close to the 
parent star, could have migrated outward (e.g. due to 
planet-planet interactions), and continue the DD burn¬ 
ing process while losing the blanketing effect from irradi¬ 
ation. These planets may appear much less inflated than 
the parked ones but should nevertheless carry with them 
signatures of DD fusion such as an eroded core. In either 
case, in our model, this implies a reduced core in today’s 
Jupiter which seems to be in agreement with recent mod¬ 
els of the interiors of today’s Jupiter (e.g. Guillot et al. 
2004; Wilson & Militzer 20I2a&b). 

From the above points, and the arguments presented 
in the previous section, we conclude that screened DD 
fusion is a viable mechanism that can provide possible 
explanations for the two key facts listed in §1. We there¬ 
fore suggest that screened DD burning operating in the 
deuterated core of giant planets is a plausible internal 
heat source and potential mechanism for the observed 
inflated radius. Based on estimates of the energy pro¬ 
duction in DD fusion at typical core temperatures, the 
amount of inflation and the trend with Tgq is broadly 
consistent with scaling relations extracted from observa¬ 
tions. Scaling relations between radius expansion and 
mass are more subtle and requires further investigation 
in our model. Seismological models of Jupiter using iner¬ 
tial modes (e.g. Dintrans & Ouyed 2001) can help probe 
the interior of the planet which hopefully would reveal 
signatures of a simmering DD burning in a reduced and 
slowly eroding core. 
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